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Abstract 

We investigate the role of two- and three-body A-nucleon forces by computing the ground state of a few A-hypernuclei 
with the Auxiliary Field Diffusion Monte Carlo algorithm. Calculations have been performed for masses up to A =41, 
including some open-shell hypernuclei. The results show that the use of a bare hyperon-nucleon force fitted on the 
available scattering data yields a consistent overestimate of the A-separation energy B A . The inclusion of a hyperon- 
nucleon-nucleon interaction systematically reduces Z?a, leading to a qualitatively good agreement with experimental 
data over the range of masses investigated. 
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1. Introduction 

The outer core of a neutron star (NS), where matter is supposed to have a density close to the saturation value, 
can be safely modeled as a gas of neutrons with a small fraction of protons and leptons (see for example Ref. [1]). 
At the larger densities reached in the inner core, the composition of matter is instead rather uncertain. The available 
observations of NS masses of order 2M Q [2, 3] has greatly reduced the softening of the equation of state (EOS) 
needed to yield a maximum mass in agreement with observational data. This fact seems to strongly question one 
of the potential sources of softening, i.e. the onset of heavier hadrons, such as hyperons, for large enough values of 
the chemical potential of NS matter. Most recent calculations on neutron matter with the inclusion of hyperons give 
a rather strong softening of the EOS [4-8]. However, results are strictly model dependent. A number of hyperon- 
nucleon interactions have been used (Nijmegen, Jiilich, ^EFT), with no overall agreement in the results. Some other 
approaches exist, that instead predict very small changes to the neutron star structure once the hyperons are included 
in the model [9, 10]. 

The only reliable benchmarks on hyperon-nucleon model forces can be made on the experiments with hyper- 
nuclei. Most of the limited available data concern single A-hypernuclei. Only a few A-nucleon scattering data are 
available. Binding energies, excitation energies and hyperon-separation energies for approximately 40 A-hypernuclei 
are available (for recent results see e.g. [1 1-13]). The lack of data reflects in the large arbitrariness in the definition of 
a hyperon-nucleon potential. 

In this paper we focus on a realistic A-nucleon potential based on the Argonne-like interaction proposed by 
Bodmer, Usmani and collaborators ([14, 15] and references therein). The potential is projected in coordinate space 
including both the AN the and AAW channel, and a realistic description of the hard-core repulsion between baryons. 
The final target will be refitting the parameters of the hyperon-nucleon-nucleon force in order to reproduce the exper- 
imental data on the A-separation energy on a set of selected hypernuclei in a rather wide mass interval. 

The A separation energy B A , defined as the difference between the binding energy of the nucleus A ~ l Z and that 
of the corresponding hypernucleus A A Z, is here computed by means of the Auxiliary Field Diffusion Monte Carlo 
(AFDMC) method [16]. 
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2. Hamiltonians and the A-nucleon interaction 



We describe nuclei and A-hypernuclei as systems made of non-relativistic point-like particles interacting via the 
following Hamiltonians: 

A-l p 2 A-l p 2 A-l A-l 

Hmc = Yj 2rL + 2 Vy ' ///,v " = H "" c + 2it + ^ VAi + ^ VAiJ ' (1) 

i=l " i<j A i=l i<j 

where A is the total number of baryons, nucleons plus the A particle. In the nucleon Hamiltonian we use the Argonne 
AV4' and AV6' two body-potentials [17], that are simplified versions of the Argonne AV18 [18]. As a test of the 
sensitivity of B A on the choice of the AW interaction, we also performed calculations with the Minnesota potential 
[19]. The A-nucleon interaction is here modeled by the two- and three-body potentials developed by Usmani et al. 
[14, 15]. It is important to note that the cited nuclear potentials do not provide the same accuracy as AV18 in fitting 
AW scattering data. In addition, three-body NNN forces are purposedly disregarded for technical reasons related to 
the AFDMC algorithm used. We rely on the assumption that when taking the difference between the binding energies 
of a nucleus and the corresponding A-hypernucleus, most of the nucleon interaction energy cancels. We shall see that 
this assumption is consistent with our results, thereby confirming that the specific choice of the nucleon Hamiltonian 
should not significantly affect the results on B A . 

Effects of the AAW force should instead not be neglected in the Hamiltonian. The lowest order in the A-nucleon 
interaction involves the exchange of two pions (TPE), in the two- as well as three-body channel (see Fig. 1). This is 
due to isospin conservation which allows only A7rE vertex and so forbids the one pion exchange process. 



AN NAN NAN AN 




AN NAN NAN AN 



(a) (b) (c) (d) 

Figure 1: Schematic construction of the A-nucleon interaction. 2n exchange processes between nucleons and the A particle. 1(a) represents the 
two-body AN channel. 1(b), 1(c) and 1(d) are the three-body ANN channels. For more details see Refs. [14, 15]. 



The AN potential is written as a sum of a central term, which also includes a A-N exchange accounting for the one 
kaon exchange process, and a spin-dependent part. The TPE process is modeled with the usual one-pion exchange 
potential acting twice, as in the Argonne forces. More details can be found in Refs. [14, 15]. 

The three-body A-nucleon interaction takes contributions from the TPE diagrams in the s-wave channel (v s 



Fig. 1(b)) and in the p-wave channel (y p 
central repulsive term v? 
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Fig. 1(c)). Short range effects are accounted for by a phenomenological 
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The definition of the functions Z(x) and Xj A as well as the set of parameters of the two-body AA^ and three-body AAW 
interactions can be found in Refs. [14, 15]. 
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We solve the ground-state of the many-body nuclear and hypernuclear Hamiltonians by means of the AFDMC 
method, originally introduced by Schmidt and Fantoni [16]. By sampling configurations of the system in coordinate- 
spin-isospin space using Monte Carlo algorithms, we evolve an initial trial wave function ¥7- in imaginary-time. In 
the t — > 00 limit, the evolved state approaches the ground-state of H. Expectation values are computed averaging over 
the sampled configurations, and for large r, the ground state energy of the system is obtained. For more details see 
e.g. Refs. [20, 21]. The extension of the AFDMC method to A-hypernuclei is technically straightforward. 

3. Results and discussion 

In Tab. 1 the results for in A He are reported. Results were obtained using different nuclear Hamiltonians. 
The second column corresponds to calculations using the AN interaction only. The numerical results with the AV4' 
and Minnesota NN potentials are partially consistent. The value obtained with AV6' has a statistically significant 
deviation. All results for Z?a are overestimated, being roughly twice the experimental value. The third column refers 
instead to calculations in which the three -body ANN force, with the set of parameters described in Ref. [14], was 
included. The three-body force significantly reduces the A-separation energy, yielding a more realistic value. It is 
very interesting to notice that the disagreement between calculations performed using different nuclear interactions is 
also reduced below statistical significance. The same conclusions were obtained by studying the l JO hypernucleus, 
for which the discrepancy between Z?a computed using different nuclear interactions is few percent. For this reason 
we are confident that the A-separation energy is not too sensitive to the details of the nuclear interaction. 



Nuclear potential 
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5.0(1) 
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6.3(1) 


5.1(1) 


3.12(2) 
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Table 1: A-separation energies (in MeV) for ^He obtained using different nuclear potentials. B^ B refers to the A-separation obtained with the 
two-body AN interaction only. B^ S+3S is the result for Ba when the three-body ANN interaction is included. In the last column the experimental 
value of B\ is from Ref. [22]. 

The A-separation energy as a function of A~ 2 ^ 3 is shown in Fig. 2. Z?a is calculated using the two-body AN 
potentials only (upper red curve) and including also the three-body ANN interaction (central blue curve). The nuclear 
potential is the AV4'. Qualitatively similar results were obtained using the AV6' or Minnesota NN forces. The dashed 
green curve is the guide line following the experimental data. It is evident that the two-body A-nucleon interaction 
alone gives an unrealistic separation energy. The saturation energy, which is estimated to be around 30 MeV, cannot 
be reproduced. The inclusion of the three-body hyperon-nucleon interaction drastically improves the theoretical 
prediction in the separation energies, and the extrapolated saturation energy becomes reasonably close to experimental 
data. We should point out that a fine tuning of the ANN parameters might further improve the AFDMC results. In 
particular, following Ref. [14] the s-wave channel has not been included in the present calculations. Preliminary 
results on A He and l JO reveal that this contribution is indeed sub-leading, and not scaling with the baryon number. 
The leading term in the ANN potential is the repulsive v A y, which is strictly dependent on A. Calculations for l JO 
performed with different values of W D , show that it is possible to move the A-separation energy closer to the expected 
result of 13.59 MeV [23], and have a better description of the A He hypernucleus. 

From these results it is clear that the role of the ANN force is important in determining the correct A separation 
energy, and should not be neglected. The fact that the leading contribution to the three-body interaction is strictly 
repulsive in the range of hypernuclei studied, and the fact that it scales with the number of baryons, suggests that the 
A-nucleon potential discussed in this paper, when applied the study of the homogeneous medium, should lead to a 
stiffer EOS for the A-neutron matter. This fact might eventually reconcile the onset of hyperons in the inner core of a 
NS with the observed masses of order 2 M . A study along this direction is in progress. 
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Figure 2: A-separation energy as a function of A -2 ' 3 . Red circles refer to the AFDMC results for the nuclear AV4' potential plus the two-body AN 
interaction. Blue diamonds are the results obtained with the same core nuclear potential but with both the two- and three- body A-nucleon forces. 
Green solid circles refer to the experimental values. 
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